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Abstract 


We investigate the stability of black holes in the viable model of f{R) = 
R+R^ gravity which was known to be the best fit for inflation. These include 
Schwarzschild and Kerr black holes. Instead of studying the fourth-order lin¬ 
earized equation around the black hole background, we use the corresponding 
tensor-scalar theory of the Starobinsky model to perform their stability. The 
Schwarzschild black hole is stable, while the Kerr black hole is unstable be¬ 
cause of superradiant instability. 
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1 Introduction 


The f{R) gravity [H El El S] has much attentions as a strong candidate 
for explaining the current and future accelerating phases EE]. On the 
other hand, f{R) black holes have included the Schwarzschild-de Sitter black 
hole [7] and Schwarzschild-anti de Sitter black hole [S] . The trace of energy- 
momentum tensor must be zero to obtain a constant curvature black hole 
when f{R) gravity couples with other matters of the Maxwell field [8], the 
Yang-Mills field [9], and a nonlinear Maxwell held |in] . 

Interestingly, f{R) = R + R^/{QM‘^) gravity [HI [121 [IS] has shown a 
strong evidence for inhation to support recent Planck data |T3]. An important 
feature of this model indicates that the inhationary dynamics were driven by 
the purely gravitational sector R^ and the scale of inhation is linked to the 
mass parameter However, it cannot work as a successful model for 

explaining late-time acceleration. 

Most of astrophysical black holes are considered to be a rotating black 
hole [H] . The stability analysis of the rotating Kerr black hole is not a routine 
work as one has performed the stability analysis of a spherically symmetric 
Schwarzschild black hole [IHl[l9lEn] because it is an axis-symmetric spinning 
black hole. The Kerr black hole has been proven to be stable against a 
massless graviton [HI [221 [23] and a massless scalar [H]. However, there 
exist another instability of the superradiance known as the black-hole bomb 
when one introduces a massive boson like scalar [251 EHl EH EH EHl EH] and 
vector 

It was hrst noted that the Kerr black hole obtained from f{R) = R + hR? 
is unstable since it could be transformed into the Brans-Dicke theory [32] . 
The Kerr solution could be obtained from a limited form of f{R) = aiR + 
a 2 R? + a^R^ + ■ ■ ■ gravity [SB]- Importantly, a perturbed Kerr black hole 
could distinguish Einstein gravity with two degrees of freedom (DOF) from 
f{R) gravity with three DOF [Mj. However, it is worth noting that the 
stability analysis of /(i?)-rotating black hole was not completely performed 
because the linearized equation for f{R) gravity contains fourth-order deriva¬ 
tive terms. One way to avoid this difficulty is to transform the limited form of 
f{R) gravity into a scalar-tensor theory with two auxiliary helds, leading to 
that the /(i?)-rotating black hole is unstable against a massive scalar pertur¬ 
bation in the Jordan frame [35] . Further, the linearized Ricci scalar equation 
obtained from the limited f{R) gravity has shown a superradiant instabil¬ 
ity if the linearized Ricci scalar is considered as a massive spin-0 graviton 
propagating on the Kerr spacetime [36] . 

In this work, we wish to focus on performing the stability of Schwarzschild 
and Kerr black holes in the specific model of f{R) = R + R?/ (6M^) gravity 


1 


because its scalar-tensor theory was clearly shown to be the Starobinsky 
model in the Einstein frame which was extensively investigated as a promising 
single held inhation model. Even though the Starobinsky potential takes the 
same form, its role in the black hole physics differs from the inhation. 

2 f{R) black holes 

We start with a specihc f{R) gravity 

with = SttG = 1/Mp. Here the mass parameter is chosen to be a 
positive value, which is consistent with the stability condition of f"{R) > 
0 [2]. The Einstein equation takes the form 

R,J'iR) - ^g^JiR) + - v^v.) f{R) = o, ( 2 ) 

where the prime (') denotes the diherentiation with respect to its argument. 
It is well-known that Eq.(|2]) provides the Kerr black hole solution with R = 
Rfiu = 0. Hereafter we denote the background quantities with the overbar. 
In this work, we use the Boyer-Lindquist coordinates to represent an axis- 
symmetric Kerr black hole with mass M and angular momentum J [TB] 

,9 2Mr\ , 5 2Mrasm'^ 9 „ , , , 

'^'^Kerr = Q^udx^dx = “ (^1--2- 2dtd(j) 

^2 , oTCr^^2 ^-,^2 a.. 

sin^ 9 d(j)^ 

(3) 

with 

A = -I- — 2Mr, cos^ 9, a = -i. (4) 

M 

In the non-rotating limit of a —)■ 0, (jS]) recovers a spherically symmetric 
Schwarzschild solution 

, 9 / 2M\ , 9 dr'^ 9 ,^9 , , 

dssch = - (}- —)dt ( 5 ) 


+^dr2 + pW+ ( 
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while the limit of a —)■ 1 corresponds to the extremal Kerr black hole. From 
the condition of A = 0{g^^ = 0), we determine two horizons which are located 
at 

r± = M ± (g) 

In the non-rotating limit ([S]), the event horizon is given by 

?"eh = 2M. (7) 

The angular velocity at the outer (event) horizon takes the form 

a a 


n = 


2Mr 


+ ' + 


ri -|- 


( 8 ) 


3 Black holes in the Starobinsky model 


Since the f{R) gravity provides three DOF, one could represent it as a scalar- 
tensor theory by introducing one auxiliary field '0 na 

Sa = j (9) 

where the superscript J means the Jordan frame. Integrating out the held 
V’ leads to the original f{R) gravity ([T]). Employing the conformal trans¬ 
formation and redehning the scalar held {ip (p) to arrive at the Einstein 
frame 

1 

—V p V 3 Mp — - 


Q-^ —y e ^ ^ 


1 + 


21/) 

MMp 


( 10 ) 


we obtain the Starobinsky model 


/ r 1 

d^x,/^[-^R - -d,4>d^ct> - V{ct>) 

with the Starobinsky potential (see Figure 1 for its graphical form) 






T <1, \ 2 
3 Mp 


The Einstein and scalar equations are given by 

1 


a,. = r* = + v) 


V' = 0, V'= \l -MpM"e 


2 - 


2 cp 

3Mp ( 1 


( 11 ) 


( 12 ) 


( 13 ) 


( 14 ) 


In the case of 0 = 0(K = 0), we obtain the Kerr solution ([3]) and Schwarzschild 
solution dS]) to ffT3)l and flTT)) . We note here that a plateau oiV |MpM^ 
appears for 0 3> 1, which was used to dehne a slow-roll inhation. 
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Figure 1: Starobinsky potential with 3MpM^/4 = 1. In the case of 0 = 
0(1^ = 0), we obtain the Kerr and Schwarzschild black hole solutions, while 
for 0 S> 1 it is sufficiently flat {V ~ 1) to ensure slow-roll conditions for 
inflation in agreement with the Planck data |14j . 


4 Linearized equations 


We start with the metric perturbation around the Kerr black hole to study 
the linear stability of the black hole 

The Taylor expansions around .R = 0 are employed to dehne the linearized 
Ricci scalar [37] as 


f{R) = f(0) + f(0)5R{h) + ---, (16) 

f{R) = f\0) + no)6R{h) + --- (17) 

with /(O) = 0, /'(O) = 1, and /"(O) = 1/3M^. The linearized equation to 
(j2|) is given by 


SR^iu{h) + 


3M2 I 




3M2 


+ V") - 


6R{h) = 0, 


where the linearized Ricci tensor and scalar are given by 

1 


SR,Ah) = 2 


6R(h) = - V^h. 


( 18 ) 


( 19 ) 

( 20 ) 


When using flT9|) and 


the linearized equation 


becomes a fourth- 


order differential equation with respect to the metric perturbation Ob¬ 
viously, it is not a solvable equation 
form 
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Another expression for flTS]) takes the 




3M2 




9fj.i' 


SR 


( 21 ) 
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whose trace equation leads to the linearized Ricci scalar equation [SB] 

(V 2 - M‘^)5R = 0 . ( 22 ) 

Choosing the Lorentz gauge of = V^h/2 and using the trace- 

reversed perturbation of = h^y — hg^y/2, Eq. ffTSjl takes a simple from [3l| 

V + ‘^R/ipuah^'" - j = 0. (23) 

Also, one could not solve fl23l) directly for 7 ^ 00 because it is still a fourth- 
order coupled equation for h^i, and h. However, its trace equation can be 
simplihed into a factorized form for h 

(v^ - M'^yh = 0 (24) 

which implies two second-order equations 

= 0, (25) 

(V 2 - M^)h = 0. (26) 

On the other hand, two linearized equations from fll3p and fll4p take the 
simple forms with = 0 and 6R = 0 

5R^,{h) = 0, (27) 

(V 2 - m2 )(^ = 0. (28) 

We note that two Eqs. fISBp and fISSp are the same but tensor equation fl2Sp 
[dH] is quite different from the linearized Einstein equation (1271) . This im¬ 
plies that the complexity of a fourth-order coupled equation oa can be 
reduced to two decoupled second-order equations fl27l) and fl28p if one em¬ 
ploys the conformal transformation and redefinition of scalar field after in¬ 
troducing the auxiliary formalism, arriving at a canonical scalar 0 with the 
Starobinsky potential in the Einstein frame. This describes a process of 
[R2 ^ R ^- 3^72 ^ - R], 

5 Stability of Schwarzschild black hole 

It is a formidable task to perform stability of the non-rotating Schwarzschild 
black hole (|B|) when one uses the fourth-order coupled equation fl25p . Actu¬ 
ally, the fourth-order derivatives appear because the tensor perturbation 
is coupled to its trace h. In addition, the Lorentz gauge condition makes the 
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stability analysis difficult because one needs to choose the Regge-Wheeler 
(RW) gauge. This amounts to double gauge-hxings. Hence, we must use 
the linearized equation flTSU to analyze the black hole stability with the RW 
gauge. As was mentioned previously, this task is not available to be per¬ 
formed because Eq. flTSU is a fourth-order differential equation with respect 
to One way to perform the stability of the Schwarzschild black hole is 
to use two Starobinsky’s linearized equations and In this case, it is 
important to note that the tensor perturbation h^i, is completely decoupled 
from the scalar ip. The stability analysis based on fl271) corresponds to that 
of the Schwarzschild black hole in Einstein gravity [THl [IHl |20]. It turned 
out that the Schwarzschild black hole is stable against the tensor perturba¬ 
tion. Furthermore, the scalar perturbation based on (1251) is stable for the 
mass-squared > 0 [57] . 

Consequently, it means that the Schwarzschild black hole is stable against 
all perturbations in the specific model of f{R) = R + R?/ 


6 Instability of Kerr black hole 

The rotating (Kerr) black hole has been proven to be stable against a mass¬ 
less spin-2 graviton [211 [221 E3] and a massless spin-0 scalar m- The sta¬ 
bility implies that normal mode solutions were allowed for tensor and scalar 
propagating on the Kerr black hole background. However, there exists an¬ 
other instability of the superradiance when one considers a massive boson 
like scalar [251 ESI EH EH ia |30] and vector [3T] around the Kerr black hole 
background. We remind the reader that either (l26l) or fl28l) is a massive scalar 
equation around the Kerr background. Here we wish to focus on the latter 
equation. Reminding the axis-symmetric background (|5]), h is convenient to 
separate the scalar into mode [55] 

r, e, 0) = , (29) 


where Sim{0) are spheroidal harmonics with —m < i < m and TZimir) sat¬ 
isfies a radial Teukolsky equation. Temporary, we may choose a positive 
frequency oj of the mode. Plugging fl29D into (1281) . one has the angular and 
radial equations for SimiO) and as 


1 

sin 6* 


de {smOdeSira) + 



M^) cos^ 9 


Adr (AdrlZemir)) = [AU - K‘^]n,m{r) 
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m 

2 n Z^-lm 

sm^ 9 


R,^ = 0, (30) 


( 31 ) 
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with U = + a^) — 2amuj + Ai^ and K = u{r‘^ + a^) — am. Here Aim 

is the separation constant whose form is given by [5^ 127] 

OO 

Aim = l{l + l) + Y.^^^^'^^^ ^ 32 ) 

k=l 


for oj M only. The Tenkolsky eqnation takes the Schrodinger form m 

- + V{r,u)i) = Eijj, i>{r) = \/~KTZ{r) (33) 

ari 

when the tortoise coordinate r* is implemented by dr* = dr and E = u‘^. 
Here, the w-dependent potential K;(r) is given by 


Vuj{r) =uj^ + 


2 , AU-K'^-M^ + a^ 


A 2 


(34) 


The approximate form of Kj — is given when keeping 1 /r-order 

2 2M(2n;2 - M^) 


E + l/^PP ^ 


r* —)■ cx) (r —)■ cx)), (35) 


while its form near the event horizon is 


— E + Vuj ^ {oj — m^lY, r* —oo (r r+). (36) 

At this stage, considering the qnalitative shape of potential K;(r*) in Fig. 
2 (see Fig.15 of Ref. [T7] and Fig.7 of Ref. [41 ] ). we could define quasibound 
states. We note that the shape of K;(r*) is slightly different from K;(r) 
because r* goes from —oo to oo while r G {r+,oo}. For this purpose, we 
impose the two boundary conditions of purely ingoing waves near the horizon 
and a decaying (bounded) solution at spatial infinity [22] . In this case, u 
should be complex because flux is absorbed at the horizon. Near the horizon 
and at the spatial inhnity, the linearized scalar takes the form m 

'ipi-o.} ~ , r* ^ -oo (37) 

V'too} = r* -)■ oo. (38) 

Since the boundary condition at the event horizon is a purely ingoing wave, 
the Kerr black holes do not admit bound states with real frequency uj. But, 
they do admit quasibound states which have complex oj = ur + ioji with 
a negative imaginary part {uj < 0), implying that the decaying held of 
jg infalling into the black holes. For the Kerr black hole, however. 
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Figure 2: Qualitative shape of the Starobinsky potential K;(r*). In the limit 
of r* —)■ oo, one hnds V{r^) —)■ a trapping potential well, a potential 
barrier, and a potential well in the ergoregion. Here, quasibound states 
appear for because the depth of potential well in the ergoregion is 

deeper than the depth of trapping potential well. 

there exists a critical frequency from flHTH : copt = 0Jc{= mVL) with up = 0, 
showing that there is no scalar flux into the black hole. For ujpi < ujc, ojj 
becomes positive, implying that the growing held is falling into black hole. 
This is the superradiant regime. It is a feature of the rotating black hole, 
but all quasibound states on the non-rotating (Schwarzschild) black hole are 
found to be decaying. Importantly, the existence of superradiant modes can 
be converted into an instability of the black hole if a mechanism to trap these 
modes in a vicinity of the black hole is provided. There are two mechanisms to 
achieve it. If one surrounds the black hole by putting a rehecting mirror, the 
wave will bounce back and forth between black hole and mirror, amplifying 
itself each time and eventually producing a non-negligible backreaction on 
the black hole background. It is not considered as a perturbation, but it 
shows a signal for instability of the black hole. Secondly, the nature may 
provide its own mirror when one introduces a massive scalar. For uj < M, 
the mass term works as a mirror effectively. 

We remember that any instability must set in via a real frequency mode 
and thus, we consider modes with |a;/| -C ur which implies ~ From 
fl55]) . a bound state of exponentially decaying mode at spatial inhnity is 
characterized by the condition 

< M^. (39) 

The three boundary conditions fl37)l - fl39|) imply a discrete set of resonances 
{un} which corresponds to bound states of the linearized scalar. 
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More precisely, according to the Hod’s argument [12] , two conditions are 
necessary to trigger the instability of the Kerr black hole when one consid¬ 
ers a massive scalar perturbation: i) The existence of an ergoregion where 
superradiant amplihcation of the waves takes place, ii) A trapping potential 
well for quasibound states should exist between the potential barrier from 
ergoregion and potential barrier from the mass (see Fig. 2). The hrst condi¬ 
tion is implemented by the superradiance condition of a; < mQ. The second 
one is supplied by the condition of the quasibound states for modes in the 
regime. This condition can be achieved by considering both fl3^ which states 
that is less than the potential hight at r = cxo and that its 

approximate derivative must be zero /dr —)■ O’*") as r —)■ cxd. Thus, one 

has the condition of 

m2 

<00^ <V^ <00^ < Ml (40) 

Combining the superradiance condition with (1101) . one hnds a restricted range 
for the mass 

M < V2u < \/2mf2 (41) 

which implies an inequality between mass M of the scalar and angular ve¬ 
locity H of the Kerr black hole 

M < V2mn (42) 

for the instability condition of the rotating black hole. On the other hand, 
the stability condition is given by 

M > (43) 


7 Discussions 

We have started with a specihc model of f{R) = R + R?/{QM'^) which is a 
fourth-order gravity theory. Considering a process of [R? —)■ R'lp — t 

— — V] have led to the Starobinsky model which is a second-order 

scalar-tensor theory in the Einstein frame. This is a famous inflation model. 
Even though the stability of Schwarzschild black hole was not carried out 
completely within the perturbed f{R) gravity, its stability was conhrmed 
within the perturbed Starobinsky model. 

In the same spirit, the stability analysis of the Kerr black hole was per¬ 
formed in the Starobinsky model because the stability analysis is a formidable 
task in the framework of the perturbed f{R) gravity. The superradiant in¬ 
stability of the Kerr spacetime is a consequence of massive modes that are 


9 


trapped inside the potential well which exists between the potential bar¬ 
rier from ergoregion and potential barrier from the mass at inhnity. In this 
case, quasibound states appear because the depth of the potential well (near 
horizon) outside the ergoregion potential barrier is deeper than the depth of 
trapping potential well. 

Finally, we wish to mention other models of fp{R) = R + \RF |13]. The 
corresponding Starobinsky potential Vp for 1 < p < 2 is steeper than p = 2 for 
large 0 (see Fig.l), while for p > 2 decreases for large 0 and it approaches 
zero [H]. Although the corresponding potentials have different shapes for 
large 0, they have the same behavior around 0 = 0 with 10,(0) = 0 which 
means that the two black holes come out as the solution. However, their 

p 1 

stability analysis seems to be unclear because 6Vp ~ pp-i -f pp-i provides 
non-integer power mass terms. The p = 2 case of our work leads to 6V' ~ p. 
In other word, we could not obtain a regular Klein-Gordon equation for fp{R) 
gravity. This is closely related to the fact that R^ for 1 < p < 2 could not 
lead to 0. 
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